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, | In this paper we find that in the thermodynamic limit and for the the ground-state normal- 

ordered ID Hubbard model the wave function of the excited energy eigenstates which span the 
Hilbert subspace where the finite-number-electron excitations are contained factorizes for all values 
of the on-site Coulombian repulsion U. This factorization results from the absence of residual energy 
interactions for the pseudofermions whose occupancy configurations describe these states. Our 
study includes the introduction of the pseudoparticle - pseudofermion unitary transformation and 
of an operator algebra for both the pseudoparticles and the pseudofermions. As the corresponding 
pseudoparticles, the cv pseudofermions (and sv pseudofermions) are 77-spin zero 2i/-holon composite 
quantum objects (and spin zero 2i/-spinon composite quantum objects) where v = 1, 2, .... The 
pseudofermion description is the most suitable for the study of the finite-energy spectral properties 
of the model. 
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I. INTRODUCTION 
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Recently there has been a renewed experimental interest in the exotic one-electron and two-electron spectral prop- 
erties of quasi-lD materials 0,0,0,3- Some of these experimental studies observed unusual finite-energy/frequency 
spectral properties, which are far from being well understood. For low values of the energy, the microscopic electronic 
(f) , properties of these materials are usually described by systems of coupled chains. For finite values of the energy 
^ ■ larger than the transfer integrals for electron hopping between the chains, the one-dimensional (ID) Hubbard model 
is expected to provide a good description of the physics of these materials 0, • This is confirmed by the recent 
quantitative studies of Refs. |l|, la, |6| . Similar unusual spectral properties observed in two-dimensional (2D) high-T c 
superconductors could result from effective quasi- ID charge and spin transport 000]. The ID Hubbard model is 
also suitable for the description of the microscopic mechanisms behind the spectral properties of the new quantum 
systems described by cold fermionic atoms on an optical lattice. Indeed, following the experimental studies of strongly 
correlated quantum systems of ultra cold bosonic atoms held in optical lattices [9j, new experiments involving cold 
fermionic atoms (such as e Li) on a ID optical lattice formed by interfering laser fields are in progress. Provided 
that the electrons are replaced by such atoms, that system can be described by the above model. However, the 
non-perturbative nature of the ID Hubbard model implies that the electronic creation and annihilation operators do 
not provide a suitable operational description for the study of the finite-energy spectral properties. Thus, the first 
step for the study of these properties is the introduction of a suitable operational description. Except in the limit of 
infinite on-site Coulombian repulsion U — > 00 ITU Il2t Il3| , the introduction of such a description is an open prob- 
£j , lem of great physical interest. For low values of energy, useful information about the effects of the non-perturbative 
electronic correlations is provided by two-component conformal-field theory [Til flit flil IrH ITsI ] . Unfortunately, that 
method does not apply for finite values of energy. 

In view of the above-mentioned finite-energy/frequency spectral properties observed in real experiments, which are 
far from being well understood, efforts towards the introduction of a suitable operational description to deal with 
the finite-ener gy p roblem are welcome. In this paper we introduce an operational representation for the ID Hubbard 
model [TgL T20L I2TL l2^| in terms of the pseudofermions. In contrast to the related pseudoparticles of Refs. [UEI, the 
pseudofermions have no residual energy interactions. We find that in the thermodynamic limit the wave function of the 
excited energy eigenstates which span the Hilbert subspace where the finite-number-electron excitations are contained 
factorizes into separated contributions corresponding to different pseudofermion branches. (A finite-number-electron 
excitation is generated by application onto the ground state of operators whose expression involves the product of a 
finite number of electronic creation and/or annihilation operators.) That factorization occurs for all values of the on- 
site Coulombian repulsion U of the ground-state normal-ordered ID Hubbard model. The pseudofermion operational 
description is closely related to the pseudoparticle representation [23L |24| , and is the natural starting point for studies 
of the finite-energy/frequency finite-number-electron spectral properties. Specifically, our starting point is a holon, 
spinon, and cO pseudoparticle representation recently introduced in Ref. 23], which refers to the whole Hilbert space 
of the model. The relation between the original electrons and these elementary quantum objects involves the concept 
of rotated electron. The rotated electrons are related to the electrons by a unitary transformation first introduced in 
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Ref. |25j . The concepts of local pseudoparticle and effective lattice widely used in this paper are introduced in Ref. 
|26| . As a result of the wave-function factorization, the pseudofermion description is more suitable for the study of the 
overlap between finite-number-electron excitations and the energy eigenstates than the pseudoparticle representation. 

As further discussed in Sec. V, the pseudofermion operational description introduced in this paper and the above- 
mentioned associated wave-function factorization are used elsewhere in the construction of a pseudofermion dynamical 
theory |27l |28| . Preliminary applications of such a theory led to a successful description of the finite-energy spec- 
tral properties observed in low-dimensional complex materials 0, 0, 01 . In spite of the absence of residual energy 
interactions, the pseudofermions are shown to be active scatterers and scattering centers in Ref. |29j| . where the 
corresponding pseudofermion scattering theory is introduced. 

The paper is organized as follows: In Sec. II the ID Hubbard model and the rotated electrons are introduced 
and useful information about the pseudoparticle representation is provided. This includes the introduction of the 
operator algebra for the pseudoparticles. The pseudofermion description and the relationship between pseudoparticle 
and pseudofermion operators are introduced and discussed in Sec. Ill, as well as the pseudofermion anticommuting 
algebra. In Sec. IV the pseudofermion energy and momentum spectra are studied and the factorization of the finite- 
number-electron Hilbert subspace of the ground-state normal-ordered ID Hubbard model is introduced and discussed. 
Finally, in Sec. V the discussion and the concluding remarks are presented. 

The pseudofermion operational description and all other issues and concepts associated with that description intro- 
duced and studied in Sees. Ill and IV correspond to a novel approach to the present quantum problem. On the other 
hand, the information about the holon, spinon, and pseudoparticle representations presented in Sec. II follows the 
results of Ref. 23] and is a review and/or different interpretation of known methods and concepts that are needed 
for the introduction of the pseudofermion description. 



II. THE ID HUBBARD MODEL, ROTATED ELECTRONS, AND THE PSEUDOPARTICLE 

DESCRIPTION 



A. THE ID HUBBARD MODEL AND ROTATED ELECTRONS 



In a chemical potential fj, and magnetic field H the ID Hubbard Hamiltonian can be written as, 



H - #SO(4) + X! ^ a ^" ' W 
a—c. s 

where the Hamiltonian H S o(a) = H H - (U/2)N + (U/4)N a has 50(4) symmetry and H H = T + U D is the 
"simple" Hubbard model. Here T = — t^2 a=1 - ^ c ] o- c i+i.ff + ^- c - i s the kinetic- energy operator and D = 

J~)j—i Cj * Cj, | Cj i Cj_ i — Yl^—i nj : i hj ! i is the electron double-occupation operator. On the right-hand side of Eq. 

/i c = 2/i, fi s = 2/io-ff, Ho is the Bohr magneton, and the number operators, S* = — \[N a — N] and S% = — ^[^t - ^jJ 
are the diagonal generators of the ry-spin and spin SU(2) algebras [soUsi! l82l | . respectively. The Hamiltonian Han u) o f 
Eq. commutes with the six generators of the 77-spin and spin SU (2) algebras and has SO (A) symmetry p30ll3ll.l32| . 
(The off-diagonal generators of these two SU(2) algebras are given in Eqs. (7) and (8) of Ref. [23j, respectively.) The 
number of lattice sites N a is considered to be large. The above electronic number operators read N = J2a=] j. N„ 
and N a = J2j=i Nj,<j, where Nj, a = c ] a c j,<? counts the number of spin-projection a electrons at real-space lattice 

site j. Moreover, the operator a (and Cj, a ) creates (and annihilates) a spin a electron at lattice site j — 1, 2, N a . 
The electronic lattice constant is denoted by a and the lattice length by L = N a a and periodic boundary conditions 
and the thermodynamic limit L —> oo are assumed. 

The momentum operator is given by P = 53 CT =t i Sfc ^(k) k where the spin-projection a momentum distribution 

operator reads N a (k) = cj, a Ck,<j and the operator c\ (and Ck,a) creates (and annihilates) a spin-projection a 

electron of momentum k. The operators cj a and Ck, a are related to the above operators cj a and cj j a by the Fourier 

transforms c\ a = [1/y/I] e lfcaj c] a and c K(y = [l/v 7 !] Ef=i ^ lkaj Cf,„ respectively. 

Throughout this paper units of Planck constant one are used and the electronic charge is denoted by — e. The 
Bethe-ansatz solvability of the ID Hubbard model JIJ is restricted to the Hilbert subspace spanned by the lowest- 
weight states (LWSs) 0, |2(J or highest-weight states (HWSs) |2lJ of the ry-spin and spin algebras, that is by the 
states whose S a and numbers are such that S a — —S^ or S a — SJ, respectively, where a = c for 77-spin and 
a = s for spin. In this paper the 77-spin and spin LWSs description of the Bethe-ansatz solution is used. In this case, 
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that solution describes energy eigenstates with electronic densities n = N/L and spin densities m = [iVj — -/VjJ/L 
in the domains < na < 1 and < ma < na, respectively. Some of the results presented in the paper correspond 
to the ranges < na < 1 and < ma < na. The description of the states corresponding to the extended domains 
< na < 1 ; 1 < na < 2 and ~na < ma < na ; — (2 — na) < ma < (2 — na), respectively is achieved by application 
onto the latter states of off-diagonal generators of the ry-spin and spin 5/7(2) algebras [23ll32T |. 

The electron - rotated-electron unitary transformation maps the electrons onto rotated electrons such that rotated- 
electron double occupation, no occupation, and spin-up and spin-down single occupation are good quantum num- 
bers for all values of U/t. The lattice occupied by rotated electrons is identical to the original electronic lat- 
tice. The electrons that occur in the ID Hubbard model are called cj a , while the operator cj a such that 
= V^(U/t) Cj a V(U/t) represents the rotated electrons, where V(U/t) denotes the electron - rotated-electron 
unitary operator. Similarly, cj a — V(U/t)cj a V^(U/t). Note that a and a are only identical in the U/t — ► oo 

limit where electron double occupation becomes a good quantum number. The unitary operators V^(U/t) and V(U/t) 
are uniquely defined for all values of U/t by Eqs. (21)-(23) of Ref. [2]|. The electron - rotated-electron unitary trans- 
formation was introduced in Ref. (25j. The rotated-electron double occupation operator D given in Eq. (20) of Ref. 
[23| commutes with the ID Hubbard model. Thus, the rotated-electron double occupation D r is a good quantum 
number for all values of U/t. 



B. THE PSEUDOPARTICLE OPERATORS 



According to the studies of Ref. |2jj, there is an infinite number of pseudoparticle branches: the cO pseudoparticles 
and the av pseudoparticles such that a = c, s and v = 1, 2, .... (Here we call cO pseudoparticles, the c pseudoparticles 
of Ref. |23| . Thus, within our notation, the general designation of av pseudoparticle refers to the av = cv branches 
such that v — 0, 1, 2, ... and av — sv branches such that v = 1,2, ....) It is shown in Ref. |2^ that for v > the cv 
pseudoparticles and sv pseudoparticles are 2^-holon and 2i/-spinon composite objects, respectively. Throughout this 
paper we follow the notation of that reference and denote the holons and spinons according to their value of 77-spin 
projection ±1/2 and spin projection ±1/2, respectively. The ±1/2 holons (and ±1/2 spinons) which are not part 
of 2i/-holon composite cv pseudoparticles (and 2^-spinon composite sv pseudoparticles) are called ±1/2 Yang holons 
(and ±1/2 HL spinons). In the designations HL spinon and Yanqholon, HL stands for Heilmann and Lieb and Yang 
refers to C. N. Yang, respectively, who are the authors of Refs. |3CM. l31| . 

The construction of the pseudofermion operational description requires the use of an operator representation for 
the av pseudoparticles, which is introduced in this section. An operational description for the pseudoparticles of 
bare-momentum q was introduced in Ref. However, that description did not take into account the holon (and 

spinon) composite character of the cv pseudoparticles (and sv pseudoparticles). Another limitation was the lack of a 
representation for the pseudoparticle operators in terms of spatial coordinates. The concepts of local av pseudoparticle 
and effective av lattice are summarized below. 

Generation and removal of pseudoparticles is in general associated with creation and/or annihilation of electrons. 
Yet there are also transitions which change the numbers of these quantum objects at constant spin-projection a 
electron numbers. Let us introduce the bare-momentum av pseudoparticle creation (and annihilation) operator w q av 
(and b q av ) which creates (and annihilates) a av pseudoparticle of bare momentum q. In addition, the local av 
pseudoparticle creation operator b\.. av and annihilation operator b Xj}0lL , are also introduced. These bare- momentum 
and local pseudoparticle operators are related as follows, 



hi - _L V p«z x i frt - _L p-"?^ ft (0) 

u q,au ~ ry / j u Xj,avi u q,au— /— / ( e u Xj,au ■ 

3=1 3=1 

The local av pseudoparticle creation (and annihilation) operator b\.. av (and b x ^ all ) creates (and annihilates) a local 
av pseudoparticle at the effective av lattice site of spatial coordinate xj = a au j, where j = 1, 2, N*„, N* u is the 
number of sites defined in Eqs. (B6)-(B8) and (Bll) and of Ref. [2j|, and 



Na L 
av au 

is the effective av lattice constant. There is an effective av pseudoparticle lattice for each av pseudoparticle branch 
|26j. This lattice has the same length L = N* v a au as the original real-space lattice. For the av pseudoparticle 
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bands such that v > 0, the discrete bare-momentum values qj are distributed symmetrically relative to zero, and are 
such that \qj\ < <7q,„. The two bare momentum values q — ±q al , limit the av pscudoparticlc Brillouin zone, where 
qav — [7t/oqi/][1 — It follows from the expressions given in Eq. (5) of Ref. [| for the around-state values of 

the number N* u that for the ground state the effective av lattice constants (01 are given by 26], 



«°o = a ; a cf = 7 j a °i = — ; a L = — > ( 4 ) 

71 •]■ TO 

where <5 = (1/a — n) is the doping concentration. The meaning of the divergences in the value of the constants aP av 
defined in Eq. J3J is that the corresponding effective av lattice has no sites, i.e. N®'* = and, therefore, does 
not exist for the ground state. This is the case of the effective cv lattices for half filling when v > and of the 
effective sv lattices for zero spin density when v > 1. ft follows that such singularities just indicate the collapse of 
the corresponding effective av lattice. This is one of the reasons why some of the expressions given in this paper refer 
to electronic densities < na < 1 and spin densities < ma < na, such that all ground-state effective av lattice 
constants Q have finite values. 

As found below, the effective av pseudoparticle and av pseudofermion lattices are identical. From use of expressions 
(@J for the ground-state effective-lattice constants a° av , one can write the ground-state number N°'* given in Eq. (5) 
of Ref. as N®'* = L/a^ v . Except for terms of order 1/L, the limiting bare-momentum values defined in Eqs. 
(B14), (B16), and (B17) of Ref. [2^| simplify and are given by q^ u = n /a°„. The ground-state expressions for these 
av pseudoparticle limiting bare-momentum values can be rewritten in terms of the original electronic lattice constant 
a as, 



9c0 = 7r / a ; Isi = k F] ; q cv = i 7r / a ~ 2 M > v>0; q° sly = [k F1 -k Fl ], v>l. (5) 

The conjugate variable of the bare-momentum qj of the av pseudoparticle branch is the above space coordinate Xj 
of the corresponding effective av lattice. This is different to the electronic operators, where the conjugate variable 
of the momentum kj is the space variable of the original electronic lattice. In reference |26|. the pseudoparticle 
site distribution configurations in the effective av lattices that describe the energy eigenstates are related to the 
corresponding rotated-electron site distribution configurations. 

The av pseudoparticle bare-momentum distribution functions N al/ (q) play an important role in the pseudoparticle 
description |23J, |24|. These functions are for all energy eigenstates the eigenvalues of the following pseudoparticle 
bare-momentum distribution operators, 



N av {q)=b\ tav b^ av . (6) 

The pseudoparticles obey a Pauli principle relative to the bare- momentum occupancy configurations, i.e. a discrete 
bare-momentum value qj can either be unoccupied or singly occupied by a pseudoparticle. Thus, the bare-momentum 
distribution functions read N au (qj) = 1 for occupied discrete bare-momentum values qj and N av (qj) — for unoc- 
cupied discrete bare-momentum values qj. Each LWS is uniquely specified by the values of the set of distribution 
functions {N ai/ (q)} such that u = 0,1, 2, ... for a — c and v = 1,2, ... for a = s. Physical quantities such as the energy, 
depend on the values of these distribution functions and numbers through the rapidity momentum functional k(q) 
and rapidity functionals A cl/ (q) and A su (q). The value of these functionals is uniquely provided by solution of the 
functional integral equations (13)-(16) of Ref. 23]. In these equations k cy (q) is the cv rapidity-momentum functional 
and the limiting bare-momentum values q^r and q av where a = c, s and v = 1,2, ... are given in Eqs. (B14), (B16), 
and (B17) of the same reference. [For the particular case of the ground state these values are provided in Eq. J5J.] 
The above integral equations correspond to a functional representation of the thermodynamic Bethe-ansatz equations 
introduced by Takahashi pp| . The rapidity-momentum functional is real and the rapidity functionals are the real part 
of Takahashi's ideal strings [20L |23|. It is useful to introduce the following cO rapidity functional A c o(q) = sinfe(g) 
where k(q) is the rapidity-momentum functional. The ground-state values of these rapidity functionals are functions 
of the bare-momentum q. Such functions are the inverse of the functions defined in Ref. |28j. 

Besides an operator representation for the pseudoparticles, another issue of importance for the introduction of the 
pseudofermion description is the definition of suitable Hilbert subspaces. For instance, an electronic ensemble space 
is a Hilbert subspace spanned by all energy eigenstates with the same values for the electronic numbers iVj and N±. 
An important concept is that of CPHS ensemble space. This is a subspace spanned by all energy eigenstates with the 
same values for the numbers {M a ,±1/2} of ±1/2 holons (a = c) and ±1/2 spinons (a — c) |23,|2J]. In general, an 
electronic ensemble space contains several CPHS ensemble spaces. Moreover, usually a CPHS ensemble space includes 
different CPHS ensemble subspaces. A CPHS ensemble subspace is spanned by all energy eigenstates with the same 
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values for the sets of numbers N c0 , {N av }, and {£ Q ,-i/2} such that a = c, s and v = 1,2,.... Here L n denotes 
the number of ±1/2 Yang holons [a = c] or ±1/2 HL spinons [a = s\. (According to the notation of Ref. [2J|, CPHS 
stands for c pseudoparticle, holon, and spinon.) 

Another tool needed for the introduction of the pseudofermion description is the ground-state normal-ordered 
pseudoparticle operational representation. As further discussed below, the finite-number-electron excitations are 
contained in a Hilbert subspace spanned by the initial ground state and excited energy eigenstates generated from it 
by processes involving changes in the occupancy configurations of a finite number of av pseudoparticles, —1/2 Yang 
holons, and —1/2 HL spinons, plus a small finite density of low-energy and small- momentum cO pseudofermion (and 
si pseudofermion) particle-hole processes for electronic densities < net < 1 (and spin densities < ma < rid). 
Throughout this paper the symbol : O : refers to the ground-state normal-ordered expression of a general operator O. 
Such a ground-state normal-ordered expression is given by that operator minus its ground-state expectation value. 
Let us then introduce the ground-state normal-ordered av pseudoparticle bare-momentum distribution operator, 

:N av (q) := N av (q) - N° u (q) , (?) 

and the —1/2 Yang holon (a = c) and —1/2 HL spinon (a = s) ground-state normal-ordered number operator, 
: L Qj _!/ 2 := 1/q,, _i/2 — L° a = L Qi -1/2- Here the operators N au {q) and L a> _i/2 are defined in Eq. © and in Eq. 
(39) of Ref. [2^, respectively, N® u (q) is the ground-state av pseudoparticle distribution function whose expression is 
given in Eqs. (C1)-(C3) of the same reference, and L° a _ 1 , 2 = is the —1/2 Yang holon (a — c) and —1/2 HL spinon 
(a = s) ground-state number. Note that the ground-state cO and si Fermi points appearing in the above-mentioned 
N® (q) and N^ 1 (q) expressions of Ref. [23| are given by, 

<7f c o = 2k F ; q° Fs i = k F i ■ (8) 

Often one can disregard the 1/L corrections and use the bare- momentum Fermi values provided in Eq. © and the 
limiting bare- momentum values given in Eq. JSJl. 

The eigenvalue equations : N au (q) : \ip) = AA r m/ (g)|^) and : £ Q ,_i/2 : \ip) — ALq, __ 1 / 2 |'0) apply to any energy 
eigenstate \ip). Here AN al ,(q) is the av pseudoparticle bare-momentum distribution function deviation and AL Qj -1/2 
is the deviation in the number of —1/2 Yang holons (a — c) or of —1/2 HL spinons (a = s). These deviations are 
given by, 

AN av {q) = N ai/ {q) - N av {q) ; AL a> _ 1/2 = i Q , _ 1/2 - i°,_ 1/a . (9) 

Such values refer to excited-energy-eigenstate deviations relative to the ground-state occupancy configurations de- 
scribed by the bare-momentum distribution functions and numbers given in Eqs. (C1)-(C3) of Ref. |25j. For these 
excited energy eigenstates, the av pseudoparticle bare-momentum distribution function and the —1/2 Yang holon 
and —1/2 HL spinon numbers read, 

N av (q)=N° u (q)+AN av (q); L a> _ 1/2 = L° a _ 1/2 + AL Q , _ 1/2 . (10) 

From use of the ground-state distribution and number values given in Eqs. (C1)-(C3) of Ref. [2j|, the operational 
relations : N cv (q) := N c „(q) for v > 0, : N s „(q) :— N sv (q) for v > 1, and : L a .-\j2 '■= i Q ,-i/2 ar e founded. Such 
relations are justified by the absence of the corresponding quantum objects in the initial ground state. 

The transition from the ground state to an excited energy eigenstate leads to a shift Q Q av jL — 0, ±n/L in the 
discrete bare-momentum value qj = [2n/L]I" 1 ' of Ref. [23j ■ That shift is such that, 

00 

OL — C. s u—1 

AN c0 + AN au odd; a = c, s , v > . (11) 

The occurrence of such a bare-momentum shift has effects on the form of the pseudoparticle operator anticommutators. 
For instance, let us consider a av pseudoparticle of bare momentum q and a a'v' pseudoparticle of bare momentum q' 
such that the bare-momentum values q and q' correspond to an excited energy eigenstate and the initial ground state, 



Q°co = E J2 AN °» even ; Q°o = ±^ 

OL — C. s u—1 

; AA^co + AN au even ; Q° = ±tt ; 



$ CtV 
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respectively. The anticommutators involving the creation and/or annihilation operators of these two pseudoparticles 
read, 



{^g, ai/i ^q' , ot'v' } — ®~ql, ol' v' $q, q' Qau/^ ^ 



QlJ 2 = :n/2, (12) 



L e+''(?-9')/2 sinQq - q'}/2) 
and the anticommutators between two av pseudoparticle creation or annihilation operators vanish. 

III. THE PSEUDOFERMION DESCRIPTION 

In this section the pseudofermion operational description and the corresponding pseudoparticle - pseudofermion 
unitary transformation are introduced. 

A. THE FUNCTIONAL CHARACTER OF THE CANONICAL MOMENTUM 

The av pseudofermion has canonical momentum q~j given by, 



L L 3 L 

Here Q% v {qj)/L is the canonical-momentum shift Junctional, 



j = l,2,...,N* v . (13) 



N* , , 

9lM = ?ij2 E *-,«v'fe,^)Aiw g/ ), (w) 

a'v' j' = l 

where the sum ^ Q , F , runs over all aV branches with finite pseudofermion occupany in the excited energy eigenstate. 
Often in this paper we use the notation av ^ cO, si branches, which refers to all av branches except the cO and 
si branches. Moreover, the summations ^2 av , J2av=c0 si' an< ^ Sai/^co si run over a ^ av branches with finite av 
pseudofermion occupancy in the corresponding state or subspace, the cO and si branches only, and all av branches with 
finite av pseudofermion occupancy in the corresponding state or subspace except the cO and si branches, respectively. 

It is shown elsewhere that the function 7r<& QI/ , a / v i (q, q') on the right-hand side of Eq. (|14|) is a elementary two- 
pseudofermion phase shift |29f. In units of 7r it is given by, 

where <& av ^ a i v > (r, r') is the corresponding rapidity two-pseudofermion phase shift expressed in terms of the variable 
r and the ground-state rapidity function hP av {q') is defined in terms of its inverse function in Ref. 28] . The rapidity 
two-pseudofermion phase shifts are the unique solutions of the integral equations i|Al(l - l|A12D of Appendix A. The 
general expression l|15f) is valid for densities < na < 1 and < ma < na. The specific two-pseudofermion phase-shift 
expressions involving the cv ^ cO branches for na = 1 (and the sv ^= si branches for ma = 0) are studied in Ref. |29| . 
The functional Q a v{<lj) / L such that, 



Q a u{q j )/2 = Qlj2 + Qt v {q j )/2, 

gives the shift in the discrete canonical-momentum value q~j that arises due to the transition from the ground state to an 
excited energy eigenstate. In turn, Q° ail {qj)/L, Eq. (|llf> . gives the corresponding shift in the discrete bare-momentum 
value qj that arises as a result of the same transition. 

The momentum shift Q® v (qj)/L of Eq. (|14|) is the part of Q av (q -)/L that refers only to the canonical momentum. It 
fully controls the following transformation which relates the av pseudofermion to the corresponding av pseudoparticle, 



(16) 
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where q~j is the discrete canonical-momentum defined in Eq. I|13fl . Note that Eq. (|13[) and the transformation 
(|16H apply both to the initial ground state and excited energy eigenstates, yet for the former state qj = q~j because 
according to Eq. (fHj) . Q% v {qj) = 0. Indeed, the pseudofermion description refers to the ground-state normal-ordered 
ID Hubbard model. Thus, there is a specific av pseudofermion description for each initial ground state. For the 
latter state the discrete bare momentum qj = [2n/L}I°" y of Eq. (Bl) of Ref. equals the discrete canonical 
momentum q~j = qj + Q^ v {qj)/L. This justifies the designation bare momentum for qj. Thus, for the ground state 
the av pseudoparticles are identical to the av pseudofermions. It follows that the ground state is invariant under the 
pseudoparticle - pseudofermion transformation and plays the role of the vacuum of the pseudofermion theory. 

The number of av pseudoparticles, N av , equals that of av pseudofermions. Moreover, we introduce the av canonical- 
momentum distribution function Af al/ (q~j) such that 

Na.v{q.j) = N a v(qj(qj)) , (17) 

where N av (qj) stands for the av pseudoparticle bare-momentum distribution function and qj — qj(q~j) is the inverse 
function of Q13[l. The function N av (qj) is the eigenvalue of the corresponding bare- momentum number operator 

A scattering theory for the pseudofermions is introduced in Ref. [2j|. In that reference it is found that the 
invariance under the electron - rotated-electron unitary transformation of the av pseudofermions corresponding to 
av pseudoparticles created at limiting bare momentum q = ±q°„ and belonging to av ^ cO, si branches implies that 
each of such cv pseudofermions (and sv pseudofermions) separates into 2v independent holons (and 2v independent 
spinons) and a cv (and sv) FP scattering center. By independent holons and spinons it is meant those which remain 
invariant under the electron - rotated-electron transformation. (The Yang holons and HL spinons are also independent 
holons and spinons, respectively.) The above designation FP stands for Fermi points. Indeed, it is found in the same 
reference that the cO and si pseudofermion and pseudofermion hole scatterers feel the created cv (and sv) FP scattering 
centers as being cO (and cO and si) pseudofermion scattering centers at the Fermi points. 

As for the pseudoparticle representation, the pseudofermion description corresponds to large values of the Hubbard 
chain length L and is thus compatible with Takahashi's string hypothesis [2(J,|23|]. However, while the pseudoparticle 
representation corresponds to the whole Hilbert space, the pseudofermion description and the associated transforma- 
tion (|16fl refer to a Hilbert subspace called pseudofermion subspace (PS). The PS is spanned by the initial ground 
state and the excited energy eigenstates generated from it by the following types of processes: 

(A) - A number N processes of finite-energy and finite-momentum pseudofermion processes such that N processes /N a — > 
as N a — ► oo, involving creation or annihilation of cO and si pseudofermions for all values of canonical momentum 
away from the Fermi points and creation of pseudofermions belonging to av ^= cO, si branches whose corresponding 
pseudoparticles have bare- momentum values obeying the inequality \q\ < q^ u . This can include a number N^ F and 
jyphNF Q £ f m }t e _ ener gy an d finite-momentum cO and si pseudofermion particle-hole processes, respectively, such that 
N% NF /N a - and N^ NF /N a - as N a - oo. 

(B) A number N processes of processes such that N processes /N a — > as N a — > oo, including creation of independent 
— 1/2 holons and independent —1/2 spinons and zero-energy and finite-momentum processes which change the number 
of cO and si pseudofermions at the corresponding Fermi points. (The latter processes transform the ground-state 
densely packed cO and si pseudofermion canonical- momentum occupancy configuration into an excited-state densely 
packed canonical-momentum occupancy configuration.) 

(C) - For densities < na < 1 and < na < ma, a number N^ F and N^ F of low-energy and small-momentum 
elementary cO and si pseudofermion particle-hole processes in the vicinity and around the canonical- momentum values 
corresponding to the bare-momentum Fermi points i~qp c0 and iqp s i of Eq. (JSJ) where t — ±1, respectively, relative to 
the excited-state densely packed av — cO, si pseudofermion canonical-momentum occupancy configurations generated 
by the above cO and si zero-energy elementary processes (B) and such that N^ lF /N a and N^ F /N a vanish or remain 
finite but small as N a — > oo. 

All PS excited energy eigenstates are generated from the initial ground state by the processes of types (A)-(C). 
Thus, the general PS excited energy eigenstates can be written as \ex) = G^(C) G'(B) G' (A)\GS), where \GS) denotes 
the initial ground state and G^(A), G^B), and G^(C) generate the processes (A), (B), and (C), respectively. Such 
generators have simple expressions in terms of pseudofermion elementary operators [271 128| . 

For simplicity, densities in the ranges < na < 1 and < ma < na are considered. The present analysis can 
be extended to other values of the densities, yet for the half-filling na = 1 or zero-magnetization ma — phases 
the excitation subspace is more reduced. The transformation associated with Eq. (|l(j[l is defined for the PS where 
it is unitary, as discussed below and in Appendix B. A crucial point of the pseudofermion theory is that the finitc- 
numbcr-clectron excitations are contained in the PS. By finite-number-electron operators, it is meant here operators 
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which can be written as a product of M electron creation and/or annihilation operators such that N /N a — ► as 
N a — > oo. The self-consistency of the pseudofermion theory is confirmed by the fact that the absolute value of the 
spectral-weight matrix elements between the initial ground state and the PS excited energy eigenstates obtained by 
use of the theory is a deceasing function of the number of finite-energy pseudofermion processes generated by the 
above operators G^(A) and G\B), vanishing as Np rnr ^ R .^ s — > oo j2^,[2g. Such matrix elements are derived by the 
pseudofermion dynamical theory reported in Refs. |27ll28| and are fully controlled by the functional Q a v{q)/2 of Eq. 
(Ilfijl . In turn, the processes (C) can involve a small finite density of elementary cO and si pseudofermion particle-hole 
processes. Indeed, the residual interactions of the corresponding cO and si pseudoparticles vanish for the subspace 
spanned by the excited states generated by these processes and thus the cO and si pseudofermion energy remains 
additive in the energy contribution of each elementary process (C) even if N^ lF /N a and N^ F /N a remain finite but 
small as N a — > oo. 

According to Eq. Ijl6j) . the processes which generate the PS excited energy eigenstates from the initial ground state 
lead to a collective canonical-momentum shift Q a v{<lj)/L — Q^ ll/ /L + Q® v (qj)/L for all the cO pseudofermions and si 
pseudofermions of the initial-state Fermi sea. In contrast to the pseudoparticles, the corresponding pseudofermions 
have no residual energy interactions. This follows from the cancelation of the av pseudoparticle residual energy 
interactions by the canonical- momentum shift functional Q^ v (q)/L of Eq. (|14|) . Such a cancelation is related to 
the form of the rapidity functionals h av (q) and rapidity-momentum functional k(q) in the PS. Introduction of the 
pseudoparticle bare-momentum distribution functions of general form given in Eq. IjlOj) in the rapidity functional 
integral equations (13)-(16) of Ref. [23| and their expansion in the small deviations ©, permits explicit solution 
of these equations. This procedure leads to expressions for the rapidity functionals A al/ (q) and rapidity- momentum 
functional k{q) in terms of the bare- momentum distribution function deviations introduced in Eq. 0. Solution of 
the above-mentioned integral equations for distributions of the general form l|l(J|) leads to first-order in the deviations 
to expressions for the rapidity-momentum functional and rapidity functionals of the following form, 



%) =k°(q(q)) ; K av {q) = k° av (q{q)) ; a = c, u = 0,1,2, ... ; a = s, ^ = 1,2, ... . (18) 

Here q(q) is the av canonical- momentum functional given in Eq. Ijl3j) with qj replaced by the continuum momentum 
q and k°(q) and h^ v (q) are the corresponding ground state functions. (These ground-state functions are fully defined 
in terms of the corresponding inverse functions in Ref. 28] . Analytical expressions for ma = and both U/t — > and 
U/t» 1 are provided in Ref. 

It is remarkable that in the PS the functionals A a „(q) and k(q) equal the corresponding ground-state functions 
hP au {q) and k°(q), respectively, with the bare momentum q replaced by the canonical-momentum functional JT3J. 
This property is behind the lack of pseudofermion residual energy interactions, as further discussed in Sec. IV. The 
canonical-momentum shift functional Ijlfijl plays a central role in the pseudofermion description of the finite-number- 
electron spectral properties. Indeed, the information recorded in the pseudoparticle interactions is transferred over to 
that functional. 

The general energy spectrum of the Hamiltonian Q depends on the quantum object occupancy configurations 
through the rapidity and rapidity-momentum functionals and can be written as, 



E = E 



SO(4) 



a— c, s 



Ha S" i E, 



SO(4) 



M c -2M c ,_ 1/2 --^ 



Na] 

2 



(19) 



where the expression of the energ y E h in terms of the rapidity-momentum functional k(q) and rapidity functionals 
A Ql /((j) is given in Eq. (20) of Ref. [2J|. As further discussed in Sec. IV, it is this functional character that is behind the 
pseudoparticle residual energy interactions. However, by re-expression of these functionals in terms of the canonical 
momentum q, the energy spectrum Ijl9jl can be written for the PS in terms of pseudofermion canonical-momentum 
distribution functions Af c o(q~j) and Af C v(q~j) as, 



N a oo N* u 

E = -2t£ JVco(fc) cosfe°fe) + 4tJ2J2^cA<jj) Re {^1 - {A° cv (qj) + ivU/4t) 2 } 



3 = 1 



u r 



M c ~Y^,2vN cv 



(20) 



The term ^ Q ^ a S" is the same as on the right-hand side of Eq. (|T§|) and N cu is the number of cv pseudofermions. 
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A crucial point of the av pseudofcrmion theory is the replacement of Eq. i|16f) of the bare-momentum q by the 
canonical-momentum q = q+Q^^q) / L. Such a procedure shows formal similarities with the usual Peierls substitution: 
The pseudofcrmion, which has no residual energy interactions, is generated from the corresponding pseudoparticle by 
the substitution of the bare-momentum q by the canonical-momentum q = q + Q® l ,(q)/L. For the PS this substitution 
renders the general energy spectrum defined by Eq. (|19|l and Eq. (20) of Ref. 24] of non-interacting form for the 
pseudofermions, as given in Eq. (|20|) . Indeed, since the bare-momentum distribution function dependent rapidity 
functional appearing in Eq. (20) of Ref. [24| are replaced by the corresponding ground-state values k°(q) and A° iy (g), 
that are independent of the set of excited-state pseudofermion canonical- momentum distribution functions {M au {q)} 
the energy l|20|) is linear in these functions. 

The form of the general energy spectrum l|20|) justifies why the shake-up effects associated with the functional 1141) 
occur in the case of the pseudofermions in the canonical momentum instead of in the energy. The dependence of 
the general energy spectrum l|20|) on that functional occurs through the canonical momentum in the argument of 
the ground-state rapidity and rapidity-momentum functions. Thus, these functions play the role of non-interacting 
spectra, since they have the same form both for the initial ground state and PS excited energy eigenstates. The shake- 
up effects associated with the two-pseudofermion phase shifts are thus felt by the pseudofermions as mere changes 
in the canonical-momentum occupancies, through the canonical-momentum shifts generated by the ground-state - 
excited-energy-eigenstate transitions. 

The pseudoparticle bare-momentum qj description is naturally provided by the Bethe-ansatz equations |23j within 
Takahashi's string hypothesis [2(j ■ We recall that the pseudoparticle discrete bare-momentum values qj are of form 
given in Eq. (Bl) of Ref. [2j| and according to Eq. (B2) of the same reference are such that qj + i — qj = 2n/L. The 
single discrete bare- momentum values qj are integer multiples of 2ir / L or of n/L |23| and bare-momentum contributions 
of order [1 / L] 3 such that j > 1 have no physical significance for the pseudoparticle description: These bare- momentum 
contributions must be considered as equaling zero. Importantly, the same is required for the pseudofermion canonical- 
momentum discrete values q~j given in Eq. I|13|) . These discrete values are also at least of the order of l/L and 
contributions of order [l/L]-' such that j > 1 must be considered as equaling zero. For instance, it is straightforward 
to find that the discrete canonical-momentum level separation, 

<U+1-QJ- T + 1 1 T . (21) 

is such that the second term on the right-hand side of Eq. I|21|) is of order [l/L] 2 , where Aq au (q) is the canonical- 
momentum shift functional given in Eq. i|14[). Thus, up to first order in l/L one finds that qj+i — q~j = 2tt/L, as for 
the corresponding discrete bare- momentum level separation given in Eq. (B2) of Ref. |23| . However, this does not 
imply that to first order in l/L the pseudofermion canonical-momentum equals the bare-momentum. Indeed, note 
that the values of the functional Q® u (qj) / L on the right-hand side of Eq. (|13fl are of order l/L and play a central 
role in the control of the finite-number-electron spectral weight distribution by the non-perturbative many-electron 
shake-up effects [27|]. What happens is that the level separation q~j + i — q~j = 2tt/L is valid locally in the discrete 
canonical-momentum space. By that we mean the following: If in the present thermodynamic limit two canonical- 
momentum values q~j and q~y differ by a small yet finite canonical-momentum value Aq = q~j — qy , then in general 
A? 7^ — j']- In contrast, for the corresponding bare-momentum values it holds that Aq = ^r\j — j']- Therefore, 
for small but non- vanishing canonical- momentum separation the difference Q® iy (qj)/L — Qavilj')/ L is not anymore 
of order [l/L] 2 and thus has physical significance. 

That only discrete canonical-momentum values qj of zero and first order in l/L are physical is an important 
property of the pseudofcrmion theory. Discrete canonical- momentum values of order [1 /L\^ with M > 1 would be 
generated by non- linear higher-order terms of the scattering phase-shift functional Q*„ (q) /2 in the pseudofermion 
canonical- momentum distribution-function deviations AA^v^g') on the right-hand side of Eq. i|ltj|) . According to 
the pseudofermion scattering studies of Ref. [23, such contributions would be associated with (Af + l)-pseudofermion 
phase shifts. However, it is shown in that reference that the corresponding (A/"+ l)-pseudofermion S matrix factorizes 
into two-pseudofermion 5* matrices. 

B. PSEUDOFERMION OPERATOR ALGEBRA 

The elementary creation and annihilation operators of the av pseudofermions can be expressed in terms of the 
corresponding operators of the av pseudoparticles as follows, 

f\j .av = Vav , av Vav j /ijj , av = Vav ^<2j . av Vav ■ (22) 
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Here V av is a unitary operator that we call the av pseudoparticle - pseudofermion unitary operator. In Appendix B 
it is shown that for the PS the operator V av which obeys Eq. (|2*2"|> is indeed unitary and given by, 

V av =CX P {^ bli^AK+Sin),^ ~ bq^, aulj] <%) =Qa„(Qj)/ L - ( 23 ) 

The canonical-momentum distribution function N av (q~j) given in Eq. I|17|) is the eigenvalue of the operator, 

■N'au (Qj ) — fqj , uv Sqj , an ■ (24) 

Keeping only the physical momentum contributions that correspond to terms up to first order in 1/L, the function 
qj (q~j ) appearing in Eq. 1|17|) is given by, 

2 N * a '"' 

Qj = ) = Qj~ A Q<™ ) = Qj ~ Y a'u> (Qj » ) (qj' ) ■ (25) 

j' = 1 

(We remind that since the functional l|16l) vanishes for the ground state, qj — q~j for that state.) On the right-hand 
side of Eq. & v a i u i(q,q') is the canonical-momentum two-pseudfermion phase shift. It is defined as, 

where q(g) is the continuum version of the function (|25() . a v (<7, <z') is given in Eq. (|15fl . ^a^.a'u'(f, r') is the 
two-pseudofermion phase shift expressed in the variable r defined by the integral equations 1A1(I - 1IA12(1 of Appendix 
A, and AP au (q') is defined in terms of its inverse function in Ref. |2q . 

Often one replaces the pseudoparticle bare-momentum summations by integrals and the corresponding discrete 
bare-momentum values qj by a continuum bare- momentum variable q. Since according to Eq. (B2) of Ref. |23| . 
the difference qj+\ — qj = 2n/L is constant for all values of j, the use of that continuum representation involves 

the replacement of Ylj=i = T^q^q al , by ^ f- q 9a " dq. In the PS, the rapidity functional A Ql/ (q) and rapidity- 
momentum functional k(q) equal the corresponding ground-state rapidity function A"„(g) and rapidity- momentum 
function k°(q), respectively, with the bare-momentum q replaced by the canonical- momentum q. It follows that in the 
PS the limiting values of the continuum canonical-momentum q are given by the ground-state limiting values iq^ v 
given in Eq. Thus, to replace the discrete canonical-momentum values by a continuum canonical-momentum 

variable q, one must replace the summations 2j=i = T^q= a - q ° by the integrals ^ f^o" dq dq ^ ■ We then introduce 
the canonical-momentum distribution function, 

X a M = '^0-M av {q)- (27) 
dq 

Here, 

« = i - y, f °"" * ^%M> . <*> 

dq at' J -ll'u' d<1 

where the function q = q(q) is given in Eq. with q~j replaced by q. The second term on the right-hand side of 
this equation is of first order in the canonical-momentum distribution function deviations. For the pseudofermion 
description only canonical-momentum and energy contributions up to first order in these deviations are physical, 
as further discussed below. As a result, for canonical-momentum distribution function deviations AJ\f a u(q) one can 
consider that, 

Atf au (q) = AN av (q), (29) 
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where in contrast to the case of Eq. I|27[) we used dq(q)/dq = 1. 
Importantly, the av pseudoparticle number operator, 



N au = ~y ] bq^avbq^av — ^ ] au^q, av = I dq 6^ Q „ 6 9i aJy , (30) 

is invariant under the pseudoparticle - pseudofcrmion transformation. It equals the av pseudofcrmion number oper- 
ator, 



^ "j fqj , av fq~j , av 
3=1 

As shown in Appendix B, this symmetry implies the unitary character of the av pseudoparticle - pseudofermion 
operator J33J). 

Moreover, the cv pseudoparticle charge and the spin and ry-spin values found and provided in Ref. ,2[| are also 
invariant under the above transformation. The same occurs for the sv pseudoparticle spin and spin-projection values 
given in that reference. The pseudoparticle - pseudofermion transformation also leaves invariant the ±1/2 Yang holons 
and ±1/2 HL spinons. The ±1/2 holon (and ±1/2 spinon) composite character of the cv ^ cO pseudoparticles (and 
sv pseudoparticles) also remains invariant under that transformation. It follows that the cv + 1 cO pseudofermions (and 
sv pseudofermions) are 77-spin zero (and spin zero) composite objects of an equal number 1/ = 1, 2, ... of —1/2 holons 
and +1/2 holons (and —1/2 spinons and +1/2 spinons). Thus, by combining Eqs. (24) and (30) of Ref. [2j| the 
±1/2 holon (a = c) and ±1/2 spinon (a = s) number operators M a ±1/2 can be written in terms of pseudofermion 
operators as follows, 



+<Za 



^ ] fq, av fq, av — ^ 



dq 



dqjq) 
dq 



fq, av fq, av ■ 



(31) 



M ai±1/2 = L at±1/2 + Y^ V ^av{q)- (32) 

v=l q=— 

Here the pseudfermion canonical- momentum distribution operator J\ au {q) is given in Eq. (|24|l and the operator 
La. ±1/2 is the ±1/2 Yang holon (a = c) and ±1/2 HL spinon (a = s) number operator provided in Eq. (39) of Ref. 
[23. Thus, for the PS all results reported in Ref. |2^| concerning pseudoparticle charge and spin transport are also 
valid for the corresponding pseudofermions. For instance, for finite values of U/t the transport of charge (and spin) is 
associated with the cO pseudofermion and cv pseudofermion quantum charge fluids (and sv pseudofermion quantum 
spin fluids). 

We recall that the bare-momentum q is the conjugate of the spatial coordinate Xj = a av j associated with the 
effective av lattice, where j = 1, 2, N*„. As for the charge (or spin) carried by the pseudoparticles and of their 
composite character in terms of chargeaons and antichargeons 23|, ±1/2 holons, or ±1/2 spinons, also the effective 
av lattice remains invariant under the av pseudoparticle - av pseudofermion unitary transformation. Indeed, the 
functional Q% v (q)/L which controls the pseudoparticle - pseudofermion transformation fTSjl does not affect the un- 
derlying effective av lattice. As shown in Ref. |29j, for the av — cO, si branches with finite occupancy in the initial 
ground state, that momentum-shift functional just imposes a twisted boundary condition. 

As for the case of the pseudoparticles, it is useful to introduce the local av pseudofermion creation operator f£. <av 

and annihilation operator f Xj ,av These operators are related to the operators f\^ av and fq.av, respectively, obtained 
from the corresponding pseudoparticle operators through the relations given in Eq. (|22|) . as follows, 



N~ 



f 1 = 
J q. olv 



V^Tv 



3 = 1 



x 3 ft 

J xj , av ' 



fq, av — 



V^iv 



Xj f 

J Xj , (XV 



(33) 



where the summations refer to the sites of the effective av lattice. The local av pseudofermion creation (and anni- 
hilation) operator /J. av (and f Xj ,av) creates (and annihilates) a av pseudofermion at the effective av lattice site of 
spatial coordinate Xj = a° av 3, where j — 1,2, ...jN*^ and aP av is the effective av lattice constant given in Eq. 
(For the PS and except for 1/L corrections we can consider that the effective av lattice constants are the ground-state 
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constants a° av .) Thus, the conjugate variable of the canonical momentum q~j of the au pseudofermion branch is the 
space coordinate Xj of the corresponding effective au lattice. The local au pseudoparticles and corresponding local 
av pseudofermions have the same effective av lattice. It follows that the local pseudoparticle and local pseudofermion 
site distribution configurations which describe the ground state and the PS excited energy eigenstates are the same. 
(These configurations are expressed in terms of rotated-electron site distribution configurations in Ref. |26j.1 

While the local av pseudoparticles and corresponding local av pseudofermions live in the same effective av lattice, 
the values of the set of discrete bare-momentum values {qf\ and canonical- momentum values {q~j} such that j 

1.2 N*^ are different and related by Eq. H25|) . There is an one-to-one relation between these two sets of discrete 

values, which keep the same order because there is no level crossing. This property follows from the values of 
the discrete bare- momentum and canonical-momentum separation given in Eq. (B2) of Ref. p3 | and Eq. I|21|l . 
respectively. 

Finally, we consider the anticommutation relations of the pseudofermion operators. It is confirmed in Refs. |2*7ll28| 
that such relations play a major role in the evaluation of finite-number-electron matrix elements between energy 
eigenstates. Let us consider the general situation when the canonical momenta q and q' of the operators /= av and 
Iq'.avi respectively, correspond to different CPHS ensemble subspaces. The anticommutator {fg av , fq',a'v'} can be 
expressed in terms of the local-pseudofermion anticommutators {/J al/ , f Xj , t a 'u' } associated with spatial coordinates 
Xj and Xj' of the effective av and a'u' lattices, respectively, as follows, 



{fq,avi fq',a'u'} — 



E 



N*, , 

a' u' 

E 

i'=i 



i(qxj-q' ay) 



{fxj,ai/i /ay , a'u'} ■ 



(34) 



The momentum operator, P = X)o-=f i Sfc -^o-(^) which is the generator for the spatial translations, com- 
mutes with the unitary electron - rotated-electron transformation. Thus, it has the same expression, P = P = 
52o-=T I Sfe N a (k)k, where N a (k) — c k a c^ a , in terms of creation and annihilation rotated-electron operators. It 
follows that the electronic lattice remains invariant under such a transformation and, therefore, the rotated-electron 
lattice and corresponding lattice constant equal those of the original electrons. Furthermore, there is a direct re- 
lation between the latter lattice and the effective au lattice populated by the local au pseudofermions associated 
with the operators /J. av and fx ,,a'v'- The above invariance implies that the local-pseudofermion anticommutators 
{fx- au> fx,, a'u'} have simple expressions. However, under the j and j' summations of Eq. (|34l) the the exotic func- 
tional character of the canonical-momentum values leads to the following unusual algebra for the au pseudofermion 
operators, 



sin ([Qauiq) - Q' au (q')]/2 



Up.V^vj-CaV^e sin([g- g -']a/2) ' W 

and the anticommutators between two au pseudofermion creation or annihilation operators vanish. Here the values 
Qavil)/ 1 ^ and Q' av (q')/2 of the canonical-momentum shift functional <|16|l refer to the CPHS ensemble subspaces 
which the canonical momenta q and q' correspond to, respectively. 

A case of particular importance is when the CPHS ensemble subspaces associated with the canonical momentum 
q' is that of the initial ground state. In that case Q' av {q')/2 — for the ground-state CPHS ensemble subspace and 
thus the anticommutation relation (|35f) simplifies to, 



i sinf Q av (q)/2 ] 

/ft f_, ,,\-S -Uq-^an iQ„uM/2 V*" KHJ ' J (oft 

XJq.auj Jq',a'u'S — u av,a'v' AT e e . . 

N* v Bin([g - q 1 ] a/2) 

Note that if sin(Q av (q)/2) would vanish the anticommutation relation l|36() would be the usual one, {/^ au , fqi^ a v'} = 
8a, a' Sv,v' Sq. q'. In contrast, in our case that quantity has in general finite values. Indeed, under nearly all ground- 
state - excited-energy-eigenstate transitions there is a canonical-momentum shift which results from non-perturbative 
shake-up effects. 

For the cu ^ cO (and su ^ si) pseudofermion branches the anticommutation relations (|35ll and i|36|) refer to 
electronic densities na < 1 (and spin densities m > 0) such that N* v /N a (and N* v /N a ) are finite. Fortunately, for 
the scattering properties only the cO and si pseudofermion anticommutation relations l|36l) are needed |27| and thus 
one can study these properties for densities such that < na < 1 and < ma < na. 
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Finally, comparison of Eqs. (|12fl and 136( 1 for the av pseudoparticle and av pseudofermion anticommutators, respec- 
tively, reveals that the absence of pseudofermion residual energy interactions implies a more complicated functional 
character for the av pseudofermion anticommutator. Indeed, for the pseudofermion description the whole information 
contained in the pseudoparticle residual interactions is transferred over to the pseudofermion anticommutator 1(36(1 . 
through the unitary transformation described by Eqs. 1(16(1 an d ((22|) . 

IV. THE PSEUDOFERMION ENERGY AND MOMENTUM SPECTRA AND THE WAVE-FUNCTION 
FACTORIZATION OF THE NORMAL-ORDERED ID HUBBARD MODEL 

In this section, we find that the description of the quantum problem in terms of the pseudofermions leads in 
the thermodynamic limit to a wave-function factorization for the PS excited energy eigenstates where the excitations 
generated by application onto the ground state of finite-number-electron operators are contained. Such a factorization 
refers to the ground-state normal-ordered ID Hubbard model. 

A. THE PSEUDOFERMION ENERGY AND MOMENTUM SPECTRA 

In the PS the energy spectrum is of the form given in Eq. (|20|l . The processes which generate the PS excited 
energy eigenstates from the ground state can be associated with two virtual excitations: (1) a finite number of ele- 
mentary processes (A) and (B) followed by small-momentum and low-energy cO pseudofermion and si pseudofermion 
particle-hole processes (C); (2) a collective canonical-momentum shift Q a ^(q)/L for all av pseudofermions and av 
pseudofermion holes with finite pseudofermion occupancy in the excited energy eigenstate. The corresponding energy 
deviation spectrum is derived from the general PS energy spectrum given in Eq. (|20|l . Such a energy deviation 
spectrum corresponds to the ground-state normal-ordered ID Hubbard model and reads, 

AE = u; + J2 A ^c0 fa ) e c0 fa ) + £ AAf sl fa )e sl fa) + ^ £ AM av fa ) e° „ fa) , (37) 

j — 1 j—l av^cO, si j — 1 

where the energy parameter luq is given by, 

lo = 2 M AM C _ 1/2 + 2/i H [AM, _ 1/2 - AJV.i] . (38) 

Here, AM a ^ 1 / 2 are the deviations in the numbers of —1/2 holons (a — c) and of —1/2 spinons (a = s) and AN s i is 
the deviation in the number of si pseudofermions. The energy deviation spectrum (|37|1 is additive in the —1/2 holon, 
— 1/2 spinons, and av pseudofermion energies. On the right-hand side of Eq. (|37H . the functions e c o(q), ^si(q), and 
e avfa are the- pseudofermion energy bands defined in Eqs. (C.15)-(C.21) of Ref. 23]. The zero-energy levels relative 
to the initial ground-state energy of these dispersions are such that, 

e c0 (±2fc F ) = e sl (±k Fl ) = &(±[tt/o - 2k*-]) = e°„(±[fc FT - fe n ]) = . (39) 

Note that the above pseudofermion energy bands equal the corresponding pseudoparticle energy bands [23L |24| pro- 
vided that the bare momentum q is replaced by the canonical momentum q. The latter bands are plotted in Figs. 6 
to 9 of Ref. Jp for ma = 0. 

In Ref. 29] it is found that for the PS both the scattering- less phase shift Q^ v /2 and the scattering phase shift 
Qavfa)/2 conserve the total energy. Thus, the virtual-excitation (2) energy spectrum vanishes and the general 
deviation-linear energy spectrum (|37l) amounts to the contributions from the excitation (1). The latter PS excitation 
involves changes in the occupancy configurations of a finite number of quantum objects generated by the elementary 
processes (A) and (B) plus the small-momentum and low-energy excitations generated by the cO and si pseudofermion 
particle-hole elementary processes (C). The excitation (2) involves a collective canonical-momentum shift of all the 
ground-state cO pseudofermions and si pseudofermions. In the thermodynamic limit, the value of the corresponding 
ground-state numbers N c o = N and 7V s i = N± approaches infinity. Thus, the self-consistency of the pseudofermion 
theory implies that the energy of the excitation (2) vanishes in that limit, so that the total energy is additive in the 
corresponding pseudofermion energies. This requirement is fulfilled because both the pseudofermion scattering events 
and the scattering-less phase shift Qq„/2 conserve the total energy [2|j. In contrast, the elementary processes (C) 
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do not involve all cO and si pseudofermions, yet the corresponding densities N^ ( /N a and L /N a can reach finite 
but small values as N a — > oo. However, we find below that in spite of these values the corresponding energy spectrum 
remains linear in the pseudofermion particle- hole numbers N^ F and F . 

The general PS energy spectrum defined by Eq. I|20l) has no pseudofermion residual-interaction energy terms 
for it is linear in the pseudofermion canonical- momentum distribution functions. In contrast, when expressed in 
terms of pseudoparticle bare-momentum distribution functions it includes residual-interaction terms. In the latter 
case it is defined by Eq. Ijl9(l and Eq. (20) of Ref. [24]]. Such residual- interaction energy terms arise from the 
dependence on the pseudoparticle bare-momentum distribution functions of the rapidity functionals k{q) and h cv (q) 
on the right-hand side of the energy expression (20) of Ref. [24j, That dependence is defined by the integral equations 
(13)-(16) of Ref. 113]. It is the surprising property that in the PS the rapidity functionals k(q) and A cl/ (q) obey 
the relations given in Eq. (|18JI . where k°(q) and A° 1/ (g) are the corresponding ground-state values, that is behind 
the absence of pseudofermion residual-interaction energy terms in expression (|20p. as discussed in Sec. III. Indeed, 
reexpression of the general energy spectrum defined by Eq. I|19|) and Eq. (20) of Ref. [24| in terms of the canonical 
momentum q leads to the replacement of the functionals k(q) and A c „(g) by their ground-state values k°(q) and 
A°„(g), respectively. This Peierls substitution procedure leads to the PS general energy spectrum l|2U|l . which is 
linear in the canonical-momentum distribution functions. The non-interacting character of that energy spectrum 
requires that in the PS the general energy-deviation spectrum must also be linear in the pseudofermion canonical- 
momentum distribution- function deviations AN au (q) of Eq. I)29[l. Indeed, if the energy-deviation spectrum involved 
non-linear contributions in the pseudofermion canonical-momentum distribution-function deviations, these quantum 
objects would have residual interactions. This implies that for the pseudofermion theory the higher-order energy 
contributions in the pseudofermion canonical-momentum distribution-function deviations have no physical meaning. 
This is consistent with full information about the spectral properties being retained in the functional Q^ u (q)/L of 
Eq. JHJ. Therefore, in order to check the consistency of the pseudofermion theory, one should confirm that in the 
thermodynamic limit only the energy-deviation spectrum (|37|) is finite, and thus that the energy terms originated by 
higher-order contributions in the pseudofermion canonical-momentum distribution-function deviations vanish in that 
limit. 

In the following we confirm that energy contributions of order equal or larger than two in the canonical- 
momentum distribution- function deviations vanish in the thermodynamic limit. For that we use the relation 
Nav{q~j) = N a v(qj(qj)) given in Eq. I|17[l. Such a relation implies that for the PS the energy contributions originated 
by processes of the same order in the canonical-momentum and bare-momentum distribution-function deviations 
differ only by the contributions of the collective shift Q®„(q)/L of the overall canonical- momentum shift Q a v(q)/L 
associated with the excitation (2). However, since the energy contributions of the excitation (2) vanish [2^, the above 
two energy contributions are identical, referring to the excitation (1) only. Thus, if the energy contributions in the 
bare- momentum distribution- function deviations of order equal or larger than two vanish in the thermodynamic limit , 
the same occurs for the corresponding contributions in the canonical- momentum distribution- function deviations. 
Moreover, we also confirm that when the ratios N^ F /N a and N^ F /N a remain finite but small as N a — » oo, the 

energy spectrum associated with the elementary processes (C) is linear in the numbers N^ L and N^ 1 F . 

For the PS both the bare-momentum and the canonical-momentum distribution-function deviations describe the 
elementary processes (A)-(C) which generate the excited energy eigenstates from the ground state. The processes 
(A) and (B) of the above excitation (1) arc in finite number. Thus, in the thermodynamic limit these processes 
involve a vanishing density of quantum objects. Let us consider a more extended Hilbert space spanned by excited 
energy eigenstates involving a small but vanishing or finite density of occupied bare-momentum values for the bare- 
momentum distribution function deviations generated by the elementary processes of type (A) and (B). In this 
case, one can derive an energy functional in terms of the bare-momentum distribution function deviations AN c o(q), 
{AN c „(q)}, and jAN s „(q)} of Eq. © where v = 1,2,.... This is achieved by solution of the integral equations 
(13)-(16) of Ref. [23 for bare-momentum distribution functions of the general form i|lU|) . Use of the obtained rapidity 
functional expressions in the general energy spectrum defined by Eq. I|19fl and Eq. (20) of Ref. |24[ leads to a general 
energy functional of the following form, 



00 




(40) 



i=l 



where the index i = 1, 2, ... refers to the order in the above bare-momentum distribution function deviations. (For the 
pseudoparticle description, the terms of order i larger than one describe the residual interactions of the pseudoparti- 
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cles.) The first and second-order terms on the right-hand side of Eq. (|4U[) are of the following general form, 

AE 1= lo q + e c0 (q)AN c0 (q) + E e sl (g) AN sl (q) + £ £ c^(?) AJV ai/ (g) , (41) 

and 

Ai? 2- Z {E E Ot(9)AJV ai/ (g) + - E ^E[Q-(^M}' (42) 

ay q— — cxu—cO. si j— ±1 

respectively. Here AJV c o(g), AN cv (q), and AN su (q) are the bare- momentum distribution function deviations given 
in Eq. © and uiq is the energy functional |J3SJ. This functional is linear in the deviations AM C .-1/2 AM s _ 1 / 2 , 
and AA s i. (The si pseudofermion number deviation of Eq. (|38|) equals the corresponding si pseudoparticle number 
deviation.) Note that the linearity in the —1/2 Yang holon and —1/2 HL spinon number deviations holds for all 
excited energy eigenstates independently of whether the values of these deviations are small or large. That behav- 
ior follows from the non-interacting character of these quantum objects j24|- The coefficients of the i — 1 linear 
terms are the av energy bands e c o(q), £si(q), and f° al ,(ci). (These bands equal the pseudofermion energy bands 
with the canonical-momentum q replaced by the bare-momentum q.) Importantly, note that the coefficients of the 
i = 2 pseudoparticle residual-interaction energy quadratic terms (14211 involve the canonical-momentum-shift functional 
Q%„(q)/L defined in Eq. (JT1J, besides the av group velocities v av (q) = de av {q)/dq and 

« Q J- For the 

pseudoparticles such a functional appears in the residual-interaction energy terms, whereas for the pseudofermions 
it is transferred over to the canonical momentum (|13(l by the unitary transformation given in Eqs. Ijlfijl an d H22J1 . 
By use of Eq. (|14|> in the first expression of Eq. I|42|l . the i — 2 energy AE 2 can be expressed in terms of suitable 
pseudoparticle / functions fau,a'v'(q,q') [13 H3 This leads to av, aV, q, and q' summations of energy terms of 
the form [1/2L] j a v,a'v' {q, q') AN au (q) AN a / y „/ (</). An important property is that the coefficients of all energy terms 
Ei of expression H4()|) of order i > 1 involve only the group velocities, the two-pseudofermion phase shifts, and/or 
bare-momentum derivatives of these functions. For instance, note that E 2 involves the velocities and the functional 
Q* l/ (q)/2 of Eq. (|14(l . This property is related to the two-pseudofermion reducibility of the (J\f + l)-pseudofermion S 
matrix [2flj . Therefore, the two-pseudofermion phase shifts control the energy expansion H4()(l for all deviation orders. 

First, we note that although the numbers N^ F and N^ F of small- momentum and low-energy elementary cO and 
si pseudofermion particle-hole processes (C) are such that the corresponding ratios N^ F /N a and N^ 1 F /N a can reach 
finite but small values as N a — > 00, such processes do not contribute to the value of the functional Q® v (qj)/2. Indeed, 
the "particle" contributions to that functional are exactly canceled by the corresponding "hole" contributions |28| . 
This behavior follows from the vanishing of the cO and si pseudoparticle residual interactions in the subspace spanned 
by the excited states generated by the elementary processes (C), as confirmed by the form of the energy E 2 , Eq. l|4*2*)) . 
and of the remaining energy terms Ei of expression l|4(J|) of order i > 2, which vanish when Q^ u (qj)/2 = 0. It follows 
that the cO and si pseudofermion particle-hole processes (C) only contribute to the leading-order energy term E\, 
Eq. (|41|l . of expression H4UI) . Thus, the energy spectrum of the small- momentum and low-energy elementary cO and 
si pseudofermion particle-hole excitations generated by the elementary processes (C) remains linear in the number 
of these processes even for a small finite density of such processes. 

Second, we address the issues related to the energy spectrum associated with the elementary processes (A) and 
(B). The energy (|40|) is an expansion in the bare-momentum distribution function deviations. For excited energy 
eigenstates whose occupancy configurations are generated by processes similar to the elementary processes (A) and (B) 
but involving a small finite density of occupied bare-momentum values for the bare-momentum distribution function 
deviations, all energy terms of deviation order i = 1,2,... are of the same order [1/L]^ 1 = L in l/L. However, 
here we are interested in the PS where the finite-number-electron excitations are contained. The PS excited-state 
deviations associated with the occupancy configurations generated by the elementary processes (A) and (B) involve 
a vanishing density of pseudofermions, —1/2 holons, and —1/2 spinons. A property of crucial importance for the 
pseudofermion theory is that for the bare-momentum distribution deviations corresponding to the PS excited energy 
eigenstates the energy terms Ei on the right-hand side of Eq. I|4(J|) associated with the excitations generated by the 
elementary processes (A) and (B) are of order [1/L] 1-1 . Therefore, for the pseudofermion-theory thermodynamic limit 
only the i = 1 energy term l|41(l is finite, all z > 1 energy terms vanishing as [1/L] 4 ^ 1 . But according to the above 
considerations, since the energy contributions in the bare-momentum distribution-function deviations of order equal 
or larger than two vanish in the thermodynamic limit, the same occurs for the corresponding contributions in the 
canonical-momentum distribution-function deviations. Therefore, the fact that for the PS occupancy configurations 
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generated by the elementary processes (A) and (B) all energy terms Ei on the right-hand side of Eq. (|40|l such 
that i > 1 vanish as [1/L] i_1 implies that for the thermodynamic-limit pseudofcrmion theory only the i — 1 leading 
energy-deviation term is finite in that limit. This is consistent with the non-interacting character of the PS general 
energy spectrum of Eqs. (121)11 for the pseudofermions. Indeed, such a non-interacting character implies that the 
energy-deviation spectrum is linear in the pseudofermion canonical-momentum distribution-function deviations and 
given by Eq. (|37(l . as discussed above. 

On the other hand, the above / functions are associated with the two-pseudoparticle residual interactions and have 
the same role as those of Fermi-liquid theory. Indeed, for small values of the energy and electronic densities n and spin 
densities m such that < na < 1 and < ma < na, respectively, the low-energy physics is controlled by the residual 
two-pseudoparticle interactions described by the i = 2 terms <|42|) of the energy functional (|4(jp. Importantly, in that 
limit the leading-order terms in 1 / L of that energy spectrum are of first order in 1/L and are contained in both the 
energy terms l|41f) and l|42[l . On the other hand, for the pseudofermions the energy spectrum (|37|) vanishes in such 
a limit, and the same information is stored in the pseudofermion canonical- momentum (|13|) through the functional 
Q% v {q)/2 of Eq. 11*4^1 . as further discussed below and in Ref. [23|. [For the pseudoparticles the functional Q% v (q)/2 
appears instead in the residual- interaction terms 1(42(1 .] 

In the low-energy Hubert subspace only the cO and si pseudoparticle branches have finite occupancies pH Il7|. 
In this case, as the limit of vanishing density of pseudoparticles contributing to the bare-momentum distribution- 
function deviations is approached, the general energy functional MUII acquires the form of the energy spectrum 
of a two-component c = cO and s = si conformal field theory [15]] . The conformal dimensions that control the 
asymptotic of the low-energy correlations functions are extracted from the finite-size energy corrections |l5j| . These 
energy corrections can be obtained by the use in the pseudo part icle energy terms 1(41(1 and 142f) of deviations AiV c o (q) 
and AN s i(q) descriptive of general low-energy excitations |l6lll7||. In spite that for the pseudofermion theory the 
energy-deviation spectrum (|3T|> has no residual-interaction terms and thus is linear in the canonical-momentum 
distribution- function deviations, it is found in Ref. |35j that the above conformal dimensions arise naturally from the 
pseudofermion momentum for low-energy excitations. 

For the pseudofermion theory, the energy l(38() on the right-hand side of Eq. 1(37(1 controls the finite-energy physics. 
The remaining energy terms refer to gapless contributions, provided that the involved pseudofermions correspond 
to the canonical- momentum values in the vicinity of the energy-band arguments of Eq. 1)39(1. For most excited 
energy eigenstates, the latter terms also lead to finite-energy contributions. The typical value of the latter energy 
contributions is of the order of the pseudofermion energy dispersion band-width per pseudofermion involved in the 
excited energy eigenstates. Note that the PS energy spectrum (13711 can have any finite value associated with the 
regions of the (k, cj)-plane where the finite-number-electron spectral functions have finite spectral weight [51 l27|. 

Provided that one considers only the contributions of first-order in the canonical-momentum distribution function 
deviations, the momentum deviation spectrum can be written in terms of the pseudofermion canonical-momentum 
distribution function deviations. It is given by, 



AP = - AM C . _ 1/a + ^ AAAcofe) ft + E E ^su%) ft + X) £ AA Mft) [- - ft] ■ (43) 
^* -i ii ii ^ 

3 = 1 v=l 3 = 1 v=l 3 = 1 

When acting onto the PS, the ground-state normal-ordered ID Hubbard model and momentum operator can be 
written in terms of pseudofermion, —1/2 Yang holon, and —1/2 HL spinon operators as follows, 



and 



H : ~ y t ^ t e ap( ( lj) '■ fqj.au '■ '■ J "qj , av 

QLU 3 = 1 



£ La,-l/2 ia,-l/2 



(44) 



N a oo 

P ■ = X ft " , cO ' : ^1 ,c0 '■ + X X ft " fli i su ' ' ' sv ' 

3 = 1 v=l 3 = 1 

oo N*„ 

+ EE^ + ^-ft] : 4- 4, -1/2, (45) 

v=l 3 = 1 a a 



respectively, where N* = N a and the operator L a _ 1 / 2 is given in Eq. (39) of Ref. 12311 . On the right-hand side of 
Eq. (|44|l . the pseudofermion energy bands are defined by Eqs. (C.15)-(C.21) of Ref. |2j|] and the —1/2 Yang holon 
and —1/2 HL spinon energies read [24| e^ c _i/ 2 = 2/i and £ls, -1/2 = 2fioH, respectively. 
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The ground-state normal-ordered Hamiltonian l|44l) and momentum operator H45|) correspond to the energy and 
momentum deviation spectra given in Eqs. (|37|l and (|43() . respectively. 

B. WAVE-FUNCTION FACTORIZATION OF THE NORMAL-ORDERED ID HUBBARD MODEL 

It is well known that both the ground state wave function and the wave function of the excited energy eigenstates of 
the ID Hubbard model can in the U It — > oo limit be constructed as a product of a spin-less fermion wave function and 
a squeezed spin wave function |lfjlll2l[l^ |. In our pseudofermion language this factorization means that in such a limit 
the expression of the momentum and energy spectra of these states is linear in the canonical-momentum distribution 
functions. It is straightforward to show that for hnite values of U/t the general energy spectrum defined by Eq. Ijl9(l 
and Eq. (20) of Ref. [24j is not linear in such functions. Therefore, the above type of factorization does not occur in 
general for the ID Hubbard model. 

Fortunately, the evaluation of finite-number-electron spectral functions can be achieved without the full factor- 
ization of the wave functions. That problem can be solved by use of the ground-state normal-ordered ID Hubbard 
model. When expressed in terms of pseudofermion operators, that Hamiltonian and associated momentum operator 
are given in Eqs. (|44|l and (|45p. respectively. In the thermodynamic limit, there is a wave-function factorization for 
the PS excited energy eigenstates. This factorization follows from the linear dependence on the canonical-momentum 
distribution-function deviations of the expressions Ij37|) and (|43|l for the energy and momentum, respectively. There- 
fore, the wave function of the energy eigenstates of the normal-ordered Hamiltonian can be expressed in the PS as 
a product of wave functions. Each wave function corresponds to a different pseudofermion branch. In excited en- 
ergy eigenstates with finite independent —1/2 holon and independent —1/2 spinon occupancy, there is also a wave 
function for these objects. In contrast, for the pseudoparticle representation the energy functional (|40l) - (|42|l includes 
bare-momentum distribution function deviation non-linear terms associated with the pseudoparticle residual interac- 
tions. For the pseudoparticle representation we cannot ignore such energy terms because they control the low-energy 
physics 0, 171. The occurrence of these energy terms mixes contributions from different branches. It follows that 
the pseudoparticle energy spectrum is not additive in the av pseudoparticle branch contributions, in contrast to 
the pseudofermion energy spectrum given in Eqs. (|37p. Thus, in this case the wave function of the excited energy 
eigenstates does not factorize in the form of a product of pseudoparticle wave functions. 

The number of wave functions contributing to the factorized wave function of a given energy eigenstate depends on 
the occupancy configurations of that state. Only the av pseudofermion branches with finite occupancy in the state 
contribute to the wave function. This contribution is in the form of a av wave function factor, as further discussed in 
Refs. [23,|28j]- The same applies to the occupancy of independent —1/2 holons and independent —1/2 spinons. 

V. CONCLUDING REMARKS 

In this paper we introduced a pseudofermion operational description for the ground-state normal-ordered ID Hub- 
bard model. We found that in the thermodynamic limit the wave function of excited energy eigenstates belonging to 
the PS where the finite-number-electron excitations are contained factorizes for all values of U/t. This factorization 
results from the absence of residual-interaction energy terms for the pseudofermions whose occupancy configurations 
describe these excited energy eigenstates. Our study included the introduction of the pseudoparticle - pseudofermion 
unitary transformation and of an operator algebra for both the pseudoparticles and pseudofermions. In the PS 
the functional Q* V (q)/L associated with that transformation exactly cancels the residual interactions of the av 
pseudoparticles, through a mechanism similar to the usual Peierls substitution. The information recorded in the 
pseudoparticle interactions is contained in that functional and is transferred over to the pseudofermion canonical 
momentum (JT3J- Moreover, we introduced creation and annihilation operators for both the pseudoparticles and 
pseudofermions and derived the anticommutation relations of these operators. Thepseudofermion anticommutation 
relations play a central role in the study of the spectral and dynamical proper ties [5|, |g, |27], |28| . 

The pseudofermion operator algebra introduced here is used in Refs. [23, |23 in the construction of a pseudofermion 
dynamical theory. That theory allows the evaluation of finite-number-electron spectral-function expressions for all 
energy values. Furthermore, the pseudofermion operational description is useful for the further understanding of the 
exotic properties displayed by low-dimensional materials. A preliminary application of the pseudofermion dynamical 
theory to the study of the one-electron spectral functions is presented in Refs. 0, 0|. The theoretical predictions 
of these references describe quantitatively for the whole finite-energy band width the one-electron removal spectral 
lines observed by photoemission experiments in the quasi-lD organic compound TTF-TCNQ. Also the predictions 
of Ref. are consistent with the phase diagram observed in a series of organic compounds. Recently, the one- 
electron problem investigated in Refs. Q, 0| was studied by the dynamical density matrix renormalization group 
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method in Ref. [3^|. The studies of the latter reference reached results similar to those of the former references. The 
pseudofermion description introduced here is also of interest for the understanding the spectral properties of the new 
quantum systems described by cold fcrmionic atoms on an optical lattice. 
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APPENDIX A: THE RAPIDITY TWO-PSEUDOFERMION PHASE SHIFTS (r,r') 

Here we provide the set of integral equations which define the rapidity two-pseudofermion phase shifts # a „ t a i v i (r, r') 
in units of 7r on the right-hand side of Eqs. (|15[) . Let us start by introducing the following Fermi surface parameters 
rj? = it sin Q/U and r® = AtB/U where the parameters Q and B are defined in Ref. j2^. In order to derive 
the integral equations which define the rapidity two-pseudofermion phase shifts Qav.a'v 1 ( r , r '), we first use in Eqs. 
(13)-(16) of Ref. |23] the bare-momentum distribution functions of the general form (|10(1 . Expansion of the obtained 
equations up to first order in the bare-momentum distribution function deviations on the right-hand side of Eqs. 
© and i|16|) leads to expression l|15|) with the two-pseudofermion phase shift ^au.a'v'ir, r') uniquely defined by the 
integral equations given below. A first group of two-pseudofermion phase shifts obey integral equations by their own. 
These equations read, 



si. cO 
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(A2) 



and 



su (r, r') = 



Here the kernel G(r, r') is given by [l7| . 
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where t(r) = — [arctan(r + rj?) — arctan(r — r°)l and l{r) = — [ln(l + (r + r° c ) 2 ) — ln(l + (r — r°) 2 )] . A second group 
of two-pseudofermion phase shifts are expressed in terms of the basic functions given in Eqs. (|A1|) - (|A3|) as follows, 
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and 
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Finally, the remaining two-pseudofcrmion phase shifts can be expressed either in terms of the functions (|A5|) - (|A7|) 
only, 

<j> c „ c0 (r> /) = I arctar Y l=lL\ I f * dr // ^co (r" r>) 

W far') = i8,,(r - r') - I T , (A9) 

2tt tt J_ r o u[l + C 1 -^-) 2 ] 

and 

^'^-;f i *' Tg^g ' (A10) 

7r7_ r o „[i + (r_i_)a] 
or both in terms of the basic functions ljAl(l - l|A3(l and of the phase shifts (|A5|) - (|A7|) . 
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and 



2?T 7T /_ r V[l + ( 2L T T-) 2 ] J- r 27r 

In the above two-pseudofermion phase shift expressions the functions 0,,, v i (x) and {x) are given in Eqs. (B5) 
and (C22) of Ref. [23j , respectively. In spite of the different notation and except for simplifications introduced here as 
a result of some integrations performed analytically, the integral equations (|Al() - ijA13l) are equivalent to the system 
of coupled integral equations (B30)-(B40) of Ref. 33] . 

APPENDIX B: THE av PSEUDOPARTICLE - PSEUDOFERMION UNITARY OPERATOR 

Here we confirm that in the PS the av pseudoparticle - pseudofermion operator V au that obeys Eq. (|22|l is unitary 
and given by expression 123|) . While that expression refers to the PS, for any energy eigenstate of the full Hilbert 
space there is a one-to-one correspondence between the set of specific discrete values of the rapidity momentum 
{kj} — {k(qj)} such that j — 1, ...,N a and the set of cO band discrete bare-momentum values {qj}. Moreover, there is 
a one-to-one correspondence between the set of specific discrete values of each cw-branch rapidity {Aj, av } = {A au (qj)} 
such that 3 = 1, ...,N* V and the set of av band discrete bare- momentum values {qj}. These branches are such that 
a = c, s and v = 1,2,3,.... This correspondence is fully defined by the integral equations (13)-(16) of Ref. [23| . 
which refer to a functional representation of the thermodynamic Bethe-ansatz equations introduced by Takahashi 
|20|. The rapidity-momentum functional is real and the rapidity functionals are the real part of Takahashi's ideal 
strings poll^. 
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Two alternative equivalent and complete descriptions for each energy eigenstate correspond to: (a) occupancy 
configurations of the above sets of discrete bare-momentum values {qj} for the av bands with finite bare- momentum 
occupancy for the state plus occupancies of the L c , -1/2 an d L s , -1/2 numbers; (b) occupancy configurations of the 
above set of discrete numbers {kj} and {Aj^ aiy } for the av branches with finite occupancy of these numbers for 
the state plus occupancies of the £ c , -1/2 an d -^s, -1/2 numbers. Since both these descriptions describe the same 
energy eigenstates, there is a uniquely defined transformation connecting the two alternative representations. That 
transformation refers to the whole Hilbert space. 

For the PS, that correspondence assumes the simple form given in Eq. I|18|) . There k°(qj) and J^ av {qj) are 
the initial ground-state rapidity- momentum and rapidity functions, respectively, whose inverse functions are given 
in Ref. [23] and q(q) is the canonical- momentum provided in Eq. I|13|) . It follows from Eq. (|18|l that for the 
PS the av pseudoparticle - av pseudofermion transformation described by Eqs. (|16[) and 122[) fully controls the 
above bare- momentum - rapidity-momentum/rapidity transformation. Since the av pseudoparticle - pseudofermion 
transformation connects two alternative representations for the complete set of orthogonality and normalized energy 
eigenstates that span the PS, the operator V av of Eq. 112211 must be unitary in the PS, as confirmed below. 

Equation (|23[1 leads to, 

Vl = exp{XX +% , )iai , - b\ j!av ] b qj , au ) ; «%■) = QtM/L ■ (Bl) 

ii 

For the PS, the functional S(qj) — Q^ v {qj)/L given in Eq. I|14fl is of the order 1/L. Thus, within the thermody- 
namic limit that the pseudofermion description corresponds to, we use the following representation for the operator 
[ &t — b\ „,.l in terms of the continuum bare-momentum o, 

[^UteW~ 6 l-]=^| 6 U- (B2) 

Next, by expanding the exponential of Eq. 1|23|) we find that to first-order in 1/L the operator (|B1(I obeys the following 
equation, 



nA^ = bl +s(qAav Vl. (B3) 



This equation is equivalent to Eq. 122fl . Therefore, this confirms that for the PS the operators (|23|) and (|B1(I indeed 
obey the relation (|22fl . Let us next show that in the PS these operators are unitary. 
The operator l)Bljl can be rewritten as, 



ftt 1 h 



(B4) 



If the operator G av is hermitian then the operator V av is unitary. Let us show that the invariance of the av 
pseudoparticle number operator l|30|) under the transformation , which implies that it equals the av pseudofermion 
number operator l|31|l . also implies the hermitian character of the operator G aL , of Eq. (|B4fl . 
Taking the transpose of the operator G av of Eq. (|B4|I leads to, 



6 L = *E & l-[ 6 



q j +S(q j ),au °qj,av\ 



(B5) 



By using the transpose of Eq. (|B2|) we rewrite the operator (|B5ll as, 



G\ v = J dq8{q)b\ <av ^-b q<au = -i^- J dq ^- 5{q)b\ aiJ 



^ +i 2^J dq dq- 



$(q) bl av b q . av 



G av - i— I dq 



(B6) 
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where 



bl 



v bq. , av] b<lj , av ■ 



(B7) 



In equations (|B6(1 and (|B6|I we have replaced the bare-momentum summations by integrals over the whole q domain 
of the corresponding av band. 

We emphasize that the operator G av of Eq. (|B4|) is hcrmitian provided that the two last terms on the right-hand 
side of the last line of Eq. ljB6(l vanish. These terms can be rewritten as, 



4 / ,,q 



r d 



dq6{q) d-q 



d r 



bq,ai> bq,ai> + i ^ 



d r 

dq dq[ 6 ^ q > av hqau 
d r " 
dq dq\^ b <i^ v K^\ + 



u g, ay 



5 r 



•A/at/ ^ay 



= o, 



(B8) 



where here the cw pseudofermion number operator l|31|) was called A7 aI , and N av is the cm/ pseudoparticle number 
operator H3()|) . Since these operators are invariant under the transformation ljl6|l they are the same operator, what 
justifies the vanishing of the operator terms l|B8(l . This shows that the operator G av of Eq. I|B4J) is hermitian and 
thus that the operator V av f23(l is unitary. 
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